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Transverse spin effects of sea quarks in unpolarized nucleons 
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We calculate the non-zero Boer-Mulders functions of sea quarks inside the proton in a meson- 
baryon fluctuation model. The results show that the transverse spin effects of sea quarks in an 
unpolarized nucleon are sizable. Using the obtained antiquark Boer-Mulders functions, we estimate 
the cos 2<f> asymmetries in the unpolarized pp and pD Drell-Yan processes at FNAL E866/NuSca 
experiments. The prediction for the cos 20 asymmetries in the unpolarized pp Drell-Yan process at 
the BNL Relativistic Heavy Ion Collider (RH1C) is also given. 
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I. INTRODUCTION 

The transverse spin phenomena appearing in high en- 
ergy scattering processes [l[ are among the most inter- 
esting issues of spin physics. Several kinds of asymme- 
tries associated with transverse spin have been observed. 
One is the single spin asymmetry in the semi-inclusive 
deeply inelastic scattering (SIDIS) processes @, H, 0, 
which requires the target nucleon to be transversely po- 
larized. Another is the cos 26 asymmetry in the unpo- 
larized Drell-Yan process [E 0] , where 6 is the angle be- 
tween the lepton plane and the hadron plane. It has 
been demonstrated that these asymmetries can be ac- 
counted for by specific leading-twist Icr-dependent dis- 
tribution functions, i.e., the Sivers function [E Q can 
explain the single spin asymmetry in SIDIS pro- 
cesses, while the Boer- Mulders function can account 
for [l2| the cos 26 asymmetry in unpolarized Drell-Yan 
processes. Despite their naive time reversal-odd prop- 
erty [H, calculations EE EI EE EE EE EE El have 
shown that these functions can be non-zero, due to the 
gauge-links [IE HH, H2| appearing in their operator defi- 
nitions. The factorization theorem [H, [24| involving In- 
dependent distribution functions has been worked out 
recently. 

Sivers or Boer-Mulders functions correspond to a cor- 
relation between the transverse spin of the nucleon or the 
quark and the transverse momentum of the quark inside 
the nucleon, respectively. Therefore, the investigation 
on them can provide information on the transverse spin 
property of the nucleon at the quark level. The experi- 
mental study of the Sivers function requires an incident 
nucleon transversely polarized, with the advantage that 
this function couples with usual unpolarized distribution 
or fragmentation functions, and therefore it can be ex- 
tracted more easily from experimental data. In the case 
of the Boer-Mulders function, it always couples with it- 
self or another chiral-odd function in the hard scattering 



Electronic 



address: 



* Corresponding author. 
mabq@phy.pku.edu.cn 

t Corresponding author. Electronic address: ivan.schmidt@usm.cl 



process. Its advantage, however, is that the spin struc- 
ture of hadrons can be studied without invoking beam or 
target polarization. In fact, the cos 26 asymmetries due 
to the Boer-Mulders function of valence quarks have been 
studied theoretically for the unpolarized pp EE HE HH 
and ttN [lEHSl Drell-Yan processes, which can be mea- 
sured in the future PAX [28| or COMPASS experiments, 
respectively. Measurements of the cos 26 asymmetries in 
the unpolarized pN Drell-Yan process can be performed 
by the FNAL E866/NuSea collaboration and at the BNL 
Relativistic Heavy Ion Collider (RHIC) [29]. This pro- 
cess is dominated by the annihilation of valence and sea 
quarks from the two incident nucleons. In this work we 
will investigate the role of the Boer-Mulders functions of 
the nucleon sea, and reveal their impact on the cos 26 
asymmetries in the unpolarized pN Drell-Yan process. 
We calculate the Boer-Mulders functions of the intrinsic 
u and d inside the proton using a meson-baryon fluc- 
tuation model. Based on the resulting sea quark Boer- 
Mulders functions, we estimate the cos 26 asymmetries 
in unpolarized pp and pD Drell-Yan processes at E866, 
and the cos 26 asymmetries in unpolarized pp Drell-Yan 
process at RHIC, respectively. The magnitude of the 
asymmetry is of several percent, and is sensitive to the 
choice of the Boer-Mulders functions of valence quarks. 
Therefore the investigation of the asymmetry in the un- 
polarized pN processes can not only provide information 
on the Boer-Mulders functions of intrinsic sea quarks, 
but also on those of valence quarks. 



II. BOER-MULDERS FUNCTIONS OF SEA 
QUARKS INSIDE THE PROTON 

The intrinsic sea quark content of the nucleon [30| is 
important for understanding its parton structure. Both 
experimental measurements and model calculations have 
shown that there is a sizable proportion of intrinsic sea 
quarks inside the nucleon, including uu, dd and ss sea 
quarks, which can not be explained by the naive quark 
model. A qualitatively successful model that accounts 
for a number of significant features of the intrinsic sea 
quarks is the meson-baryon fluctuation model (3l| , which 
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is based on the light-cone Fock state expansion of the 
nucleon. The main assumption of this model is that the 
intrinsic sea quarks are multi-connected to the valence 
quarks and can exist over a relatively long lifetime within 
the nucleon bound state, therefore the intrinsic qq pairs 
can arrange themselves together with the valence quarks 
of the target nucleon into the most energetically-favored 
meson-baryon fluctuations. It is easy to understand that 
the most important fluctuations are most likely to be 
those closest to the energy shell and thus have mini- 
mal invariant mass. Such fluctuations are necessary part 
of any quantum-mechanical description of the hadronic 
bound state in QCD and have also been incorporated into 
the cloudy bag model and Skyrme solution to chiral the- 
ories. Therefore, in the meson-baryon fluctuation model, 
the nucleon can fluctuate into an intermediate two-body 
system of a baryon and a meson which are in turn com- 
posite systems of quarks and gluons. For example, the 
proton can fluctuate into a mt + or A ++ 7r~ configuration, 
and the d or u inside tt + or ir~ can be viewed as d or u 
distributed in the proton. 

Based on the meson-baryon fluctuation model one can 
also calculate the momentum distribution of the intrinsic 
sea quarks inside the proton, from the following convolu- 
tion form 




where f M j BM {v) lii probability of finding the meson 

M in the BM state with light-cone momentum fraction y, 
and qM {x/y) is the probability of finding a quark q in the 
meson M with light-cone momentum fraction x/y. The 
function f i (y) can be determined from the form 

^m/baM = j d ' kT \*BM{y, k T )| 2 , (2) 

where ^ bm{x, ky) is the light-cone two-body wave func- 
tion of the baryon-meson system, for which we will choose 
a Gaussian or a power-law behavior, 

^bm{xM) = A Z5 cxp(-M 2 /8a| ) ), (3) 
* BU (xM = A' D (l+M 2 /a 2 D )- p , (4) 

respectively. Here M 2 = (M| f + kf,)/x+ {M B + k|)/(l - 
x) is the invariant mass squared of the baryon-meson 
system, ao sets the characteristic internal momentum 
scale, and An or A' D is the normalization constant and 
can be fixed by 

J dxd 2 \t T \^ BM {xM)\ 2 = 1. (5) 

We point out that Eqs. © and (f?]) are boost-invariant 
light-cone wavefunctions which emphasize multi-parton 
configurations of minimal invariant mass. 

Eq. JT]) can be extended to calculate the Boer-Mulders 
functions of the intrinsic u and d antiquarks inside the 



proton (denoted as h 1 ' u and h 1 ' ). In fact, accord- 
ing to the meson-baryon fluctuation model, there are 
pion components in the intermediate state of the pro- 
ton. On the other hand, as shown in Ref. 19], non- 
zero Boer-Mulders functions of valence quarks inside the 
pion (denoted by h^. q ) can been calculated using the 
quark-spectator-antiquark model. Therefore, a convolu- 
tion form similar to Eq. ([T]) can be applied to calculate 

hi and h^' d , as follows: 

The transverse momentum of the antiquark in the nu- 
cleon should be a superposition of the transverse momen- 
tum of the antiquark in the pion and of the transverse 
momentum of the pion in the nucleon. For simplicity we 
neglect here the transverse momentum of the pion. As- 
suming that the probabilities of the proton fluctuating to 
nir + and A ++ 7r~ (denoted as V p ^ n „+ and V p _,a++tt- j 
respectively) are both 12%, and using the functions 
given in Ref. p^ ]. we obtain the numerical results for 

h x ' u {x, kj-) and h x ' (x, fcf*), shown in Fig. [T] Two sets 
of functions are given, corresponding to the choice of 
Gaussian (upper panel) and power-law type (lower panel) 
wavefunctions of the meson-baryon system, respectively. 
The parameters used in the calculation are: ocd = 0.33 
GeV, P — 3.5 (as pointed out in Ref. [32], perturba- 
tive QCD predicts a nominal power-law fall off at large 
corresponding to P = 3.5), following the choice in 
Ref. [3l|. The quark mass used in the calculation is 
m = 0.3 GeV. For the consistency here we use the same 
value for the mass of the quark inside the proton and the 
pion, which is different from that in Ref. [19] where we 
use 0.1 GeV for the quark mass of the pion and 0.3 GeV 
for the quark mass of the proton. This difference on the 
value of quark mass in the pion will bring a quantita- 
tive difference on the numerical result, but the result of 
the asymmetry and the main conclusion are qualitatively 
not changed. It is interesting to point out that the Boer- 
Mulders function of u or d, calculated from the Gaussian 
type wavefunction, qualitatively agrees with that from 
the power-law type wave function. We should notice that 
this result is based on the assumption that the numbers 
of the sea quark pairs uu and dd inside the nucleon are 
the same, since we assume V p ^ n7T + = V p ^a++it- here. 
As we know that there should be an excess of dd pair 
over uu pair in the proton, which is a consequence of Got- 
tfried sum rule [331 ] , in next section we will apply another 
set of V p ^ n7T + and V p ^,a++it- j which is constrained by 
the existing parametrization that can produce the flavor 
asymmetry of the sea quarks, and investigate the conse- 
quent asymmetries in the unpolarized pp and pD Drell- 
Yan processes. Boer-Mulders functions for sea quarks 
may be obtained from other models that also allow sea 
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FIG. 1: The Boer-Mulders functions of u quark h-y '"(a;, kjj (left column) and d quark h x ' (x, kj_) (right column) inside the 
proton as two-dimensional densities. The upper panel and lower panel correspond to the choice of gaussian type and power-low 
type light-cone wave functions for the meson-baryon system, respectively. 



quarks in the nucleon, such as the chiral quark and me- 
son cloud models. In these models there are also 7r meson 
components in the Fock state expansion of the nucleon. 



III. THE cos 20 ASYMMETRIES IN THE 
UNPOLARIZED pp AND pD DRELL-YAN 
PROCESSES 



In this section, we will calculate the cos 2<p asymmetries 
in the unpolarized pp and pD Drell-Yan processes, based 
on the Boer-Mulders functions of sea quarks obtained 
above. The general form of the angular differential cross 
section for the unpolarized Drell-Yan process is 



1 da 



1 



4tt A 



(l + Xcos 2 6 + /isin26'cos0 



-— sin 6cos2(j)J 



(8) 



polarized pp Drell-Yan process is 



da(pp -> UX) 
dfldxidx2d 2 q±_ 



a 



j2O2 (l + cos 2 0) ]T <%F\fi{x uT >\) 

q—u 7 d 

Xf!(x 2 ,k 2 ± )} + (q^q), (9) 
where we have used the notation [34j |: 

F[ ■■}= J d 2 p ± d 2 k ± 6 2 (p ± + k x - q ± ) x {• • •}, (10) 

and the transverse momentum dependence of the cross- 
section is implicit. 

The corresponding cross section for the unpolarized pD 
Drell-Yan process is then 



da(pD UX) 
d£ldx\dx2d 2 q± 



+e 2 J?(x 1 ,pl))x(K(x 2 ,ki) 
+fi(x 3 ,** ± ))] + (q~q), (11) 



The <$t independent differential cross-section for the un- /" = /" + ff 



Where we have used the isospin relation 

/r, 



ID 



n 



u/p 
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Eqs. (J9j) and (jTTj) are lowest order parton model ex- 
pressions, where it is extended in order to include trans- 
verse parton momenta. 

The cos 2<j) dependent cross-sections for unpolarized pp 
and pD Drell-Yan processes, in terms of Boer-Mulders 
functions, are 



da(pp -> UX) 



COS 2<p 



d£ldxidx2d 2 c\j 

^sin 2 0cos20 Yl 4^(P^ k ^'Vi>Pi) 

q— u,d 

xhi'%x 2 ,k 2 ± )] + (q^q); (12) 
da{pD -> UX) 

cos 20 



dfldxidx2<PqA 



or 



■ sin 2 0cos2</)J r [x(pj_,k^)(e 2 1 /i] L,tl (xi,p 



12Q : 



+e > d hf" l (xi, pl))(h±*(x 2 , ki) + ^"(aa, ki)] 

+ (9 (13) 



respectively, where 

X(p±, k x ) - (2h ■ P± h • k ± - p_L ■ k ± )/Af 2 , 



(14) 



here h = qj-/ and M p is the mass of the proton. 

Therefore we can express the cos 2<f> asymmetry coeffi- 
cient defined in Eq. © as (A = 1, \i = 0) 



2E 



(15) 



2T[ X {e\h^ + e*hi' d )(hi* + ftjj] + (g ^ g) 



(16) 

where we have omitted the arguments of the distribution 
functions. 

Eqs. (fT5|) and (fT6|) give the explanation for the cos 2<p 
asymmetry observed in the unpolarized Drell-Yan pro- 
cess from the view of Boer-Mulders function. We note 
that there are other theoretical approaches that have 
been proposed to interpret this asymme try such as high- 
twist [35[ Hi| and QCD vacuum effects [33] • In Ref. [38j 
detailed comparison between the QCD Vacuum effect 
and Boer-Mulcrs effect has been done. 

The E866 Collaboration at FNAL employs a 800 
GeV/c proton beam colliding on protons or deuterons 
fixed targets to measure the unpolarized Drell-Yan pro- 
cesses pp — > /i + /i~A and pD — > /i + fi~ X [39j. The main 
goal of E866 is to study the asymmetry of the nucleon 
sea distribution [4(| El[ . These experiments can also be 
applied to study the cos 20 asymmetry of the lepton pair, 
especially the role of the sea quarks contribution to the 
asymmetry. In the following we will estimate the asym- 
metry v p and vd for the E866 experiment. In order to 
calculate the asymmetry we also need the specific form of 



Boer-Mulders functions of valence quarks inside the nu- 
cleon. Since there are no experimental measurements on 
those functions yet, we use some model results for them, 
such as those given in Ref. [l7| , within a quark and spec- 
tator diquark model of the nucleon. We will take two 
different extreme options as follows. 

• Option I: Consider the limit with only the spectator 
scalar diquark included, which means that h^' u — 



x,s 



and 



0. Here we use hf' S and h{ 



to represent the contributions to the Boer-Mulders 
functions with only spectator scalar diquark and 
vector diquark included, respectively. This option 
was also applied in Ref. [l9[ 

• Option II: Consider the case with both the specta- 
tor scalar diquark and vector diquark included. In 
this option h 1 = f % +f"i and — h l , 
and they have opposite sign. 

The scalar and vector components of the proton are [l7j 

h Y (x,k T ) = -a s N s (l-x) ——^—.(IT) 

, xv/ .as 4 Ar . ,o M (2xM + m). „. 
*t' V M) = - 3 a s N v (l- x f [L 2^2 +k 2 )3 J p ) 

respectively, where N s / V is a normalization constant, and 
L 2 S/V = (l-x)A 2 + x M s % -x(l-x) M 2 . (19) 

Here A is a cutoff appearing in the nucleon-quark-diquark 
vertex and Md is the mass of the scalar diquark. 
With the following kinematical cuts 

4.5 GeV < Q < 9GeV and Q>llGeV, 
0.1<xi<l, 0.015 < x 2 < 0.4, 

we calculate the Qr-dependent cos 2(f> asymmetries in the 
unpolarized pp and pD Drell-Yan processes at E866, as 
shown in Fig. [5] In the calculation we adopt the anti- 
quark Boer-Mulders functions hf'^ calculated from the 
Gaussian type wavefunction of the baryon-meson sys- 
tem. One observation from Fig. [5] is that the asymme- 
try calculated according to option I is of several percent, 
while that according to option II is around one percent, 
indicating that the size of the asymmetry is sensitive 
to the choice of the Boer-Mulders functions of valence 
quarks. Nevertheless, our calculation predicts a quite 
smaller cos 2<fi asymmetry in the unpolarized pN Drell- 
Yan process compared to that in the unpolarized ttN 
Drell-Yan case [g, [ijj. Again, comparing the asymme- 
tries calculated from the two different options, we find 
that Up and ud have similar sizes in option I, while those 
of v p and vd in option II can be very different, and even 
have opposite signs. Therefore, the investigation of the 
cos 2(f> asymmetry in the unpolarized pN Drell-Yan pro- 
cess, can not only give information on the Boer-Mulders 
functions of sea quarks, but also about those of valence 
quarks. 
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FIG. 2: The cos 2cj> asymmetries v v (thick lines) and Uo (thin 
lines) for the E866 experiment, as functions of Qt, with 
kinematical cuts 4.5 GeV < Q < 9GeV and Q > HGeV, 
0.1 < xi < 1, 0.015 < x 2 < 0.4. The dashed and solid lines 
show the asymmetries calculated from options I and II, re- 
spectively. 
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FIG. 3: The flavor asymmetry of the sea quark distribu- 
tions x(d(x) — u(x)). The solid curve shows the result from 
CTEQ6L parametrization at Q — 0.6 GeV. The dashed and 
dotted line correspond to the results from Set I and II, re- 
spectively. 



The main success of the baryon-meson fluctuation 
model is that it can qualitatively describe a variety of 
asymmetries of the nucleon sea distributions, such as 
the flavor asymmetry d(x) — u(x), the strange sea asym- 
metry s(x) — s(x), and the polarized sea asymmetry 
As(x) — As(x). In the previous section we have assumed 
that the probabilities of the proton fluctuating to nn + 



FIG. 4: Same as Fig. [2] but using Set II for the sea quark 
distributions which produce the flavor asymmetry of sea quark 
indicated by the Gottfried sum rule. 



and A ++ 7r~ are 12% (here we denote it as Set I for the 
sea quark distributions). From this set we calculate the 
flavor asymmetry x{d — u), as shown by the dashed line 
in Fig. [31 and compare it with that from the CTEQ6L 
parametrization at Q=0.6 GeV (shown by the solid line 
in the same figure). None-zero asymmetry of the sea 
quark densities is found in this case, but is quantitatively 
much lower than the known parametrization. In a fur- 
ther consideration, we will take into account the flavor 
asymmetry of the sea quark and adopt V p ^ n7T + = 15% 
and V p -,a++tt- = 1% (denoted as Set II), motivated by 
the expectation that there should be an excess of uu pair 
over dd pair, which is a consequence of Gottfried sum 
rule [HI]. This set has also been adopted in Ref. 
We show the asymmetry x(d — u) from Set II by the 
dotted line in Fig. [3j and find that it agrees with the 
flavor asymmetry from CTEQ6L parametrization fairly 
well. From the sea quark distributions of Set II, we give 
the cos 2(f> asymmetries in the unpolarized pp (thick line) 
and pD (thin line) Drell-Yan process at E866, as shown 
in Fig. |U Two options of the valence quark distribution 
are applied in this calculation. Again, the asymmetries 
Up and ujj have similar sizes in option I, and those in 
option II are very different. In Fig. 2J a larger negative 
v p from option II is predicted compared to that given in 
Figd 

As we don't know the exact form of Boer-Mulders func- 
tions of valence quarks, we admit that our prediction, in 
a certain aspect, relies on functions given in [l7j . where 
the Boer-Mulders functions for the u and d quarks are of 
the opposite sign. Actually different models or theoreti- 
cal considerations predict different flavor dependence of 
valence Boer-Mulders functions. A calculation based on 
the MIT bag model gave [HI h^' d — ^hf' u , while the ar- 
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to calculate the asymmetries at E866, we estimate the 
cos 2(f> asymmetry v v at RHIC with the kinematical con- 
straints y/s = 200 GeV, — 1 < y < 2, here y is the rapidity 
defined as y — glnfV The Qy-dependent asymmetries 
for Q — 4 (thin line) and Q = 20 GeV (thick line) are 
shown in Fig. [5] Here we do not average the asymmetry 
over a range of Q, as we have done for the asymmetry 
at E866, since we would like to choose two different fixed 
values of Q to study the Q~dependence of the asymmetry 
at RHIC. Since the evolution equations of unintegrated 
parton distributions, especially that of the Boer-Mulders 
functions, are unknown, we assume that they scale with 
Q. Therefore the different asymmetries at different Q 
values, as shown in Fig. is not a consequence of evolu- 
tion, but just a kinematical effect. 



FIG. 5: The Qr-dependent cos2<^> asymmetries v v at RHIC, 
with the kinematical conditions y/s — 200 GeV, — 1 < y < 2. 
The thin and thick lines show the asymmetry at Q = 4 and 
20 GeV, respectively. The dashed lines show the asymmetry 
from option I and Set I or II (in option one the asymme- 
tries in Set I and II are the same), and solid lines show the 
asymmetries calculated from options II and Set I, dotted lines 
show the asymmetries calculated from options II and Set II, 
respectively. 



gument based on large- N c showed that hf' d = hj~' u (43j 
modulo XJN C corrections. Recent lattice work [44( and 
works 45] on impact parameter representation of gener- 
alized parton distributions (GPDs) state that h 1 ' u and 
h^' are of the same sign. In those cases the asymmetries 
may be different from the prediction given in this work. 
In our recent paper [27| . it has been suggested to study 
the flavor dependence of valence Boer-Mulders functions, 
through np and nD Drell-Yan processes that can be per- 
formed in the coming years by COMPASS Collaboration. 
In this paper, We have taken two options for the valence 
Boer-Mulders functions to study how the asymmetry de- 
pends on them. Nevertheless, the measurement of the 
cos 2<p asymmetries in the unpolarized pp and pD pro- 
cesses may provide useful constraint on the Boer-Mulders 
function of valence quarks. 

Another promising test ground of cos 2<j) asymme- 
tries in the unpolarized pN Drell-Yan process is BNL 
RHIC [2{|, where the Drell-Yan process pp — > l + l~X can 
be measured. There are already numerous studies [13, HE 
H3, [48| on the single spin and double spin asymmetries 
in Drell-Yan processes at RHIC, which require at least 
one incident proton to be polarized. The unpolarized 
experiment is also feasible at RHIC. Taking the same 
sets of functions (Set I and II for sea quark distributions 
and Option I and II for valence quark distributions) used 



IV. CONCLUSION 

In summary, we have applied a meson-baryon fluctu- 
ation model to calculate the Boer-Mulders functions of 
u and d inside the proton, and provided a first estimate 
of the transverse spin effect of sea quarks in the unpo- 
larized nucleon. In the calculation we adopted both the 
Gaussian and the power-law type wavefunctions for the 
meson-baryon system. We found the sizes of h^' M or h^' d 
from these two types of wavefunctions to be qualitatively 
similar. From the obtained antiquark Boer-Mulders func- 
tions, we analyzed the cos 2<f> asymmetries in the unpo- 
larized pp and pD Drell-Yan processes at E866. Two sets 
of the sea quark distributions are applied in the calcula- 
tion, in which Set I is flavor symmetric and Set II reflects 
the sea quark flavor asymmetry indicated by Gottfried 
sum rule. The size of the asymmetries is at most sev- 
eral percent, in the case that the Boer-Mulders functions 
of valence quarks are chosen with only spectator scalar- 
diquark contributed. The size of the asymmetries, as 
well as the relative size between the asymmetries in the 
unpolarized pp and pD Drell-Yan processes, relies signifi- 
cantly on the choice of Boer-Mulders functions of valence 
quarks. We also estimated the cos 2(f) asymmetry in the 
unpolarized pp Drell-Yan process at RHIC. This investi- 
gation suggests that unpolarized pN Drell-Yan processes 
on hadron colliders are very helpful for a better under- 
standing about the role of transverse spin of both valence 
and sea quarks in the unpolarized nucleon. 
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